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Abstract: 

Photoproduction of 2-jet events is studied with the HI detector at HERA. Parton 
cross sections are extracted from the data by an unfolding method using leading order 
parton-jet correlations of a QCD generator. The gluon distribution in the photon is 
derived in the fractional momentum range 0.04 < x 7 < 1 at the average factorization 
scale 75 GeV 2 . 
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1 Introduction 



The interaction of electrons and protons at the HERA collider is dominated by photoproduc- 
tion processes: electrons scatter through small angles and emit quasi-real photons, which then 
interact with the protons. The center of mass (CMS) energies in the 7p system reach up to 
300 GeV. A fraction of these events has large transverse energy in the final state including the 
formation of jets, as has been reported in several recent publications ^. B|. 

The jet production can be well described in the framework of QCD. In this picture the 
photon couples either directly to a parton of the proton, or indirectly via the photon's own 
parton content. The first are called direct processes, which include the QCD-Compton (Fig. la) 
and photon-gluon fusion diagrams, while the latter are usually refered to as resolved processes. 
An example is shown in Fig. lb. 




Figure 1: Examples of diagrams for direct (a) and resolved photon (b) processes in electron- 
proton scattering. 

Predictions for the jet cross section are usually obtained in leading order (LO) QCD by 
convoluting the parton densities in the photon and in the proton with hard partonic scattering 
cross sections calculated at the tree level. The partons leaving the hard scattering reaction are 
identified with hadronic jets. Due to the reduced center of mass energy of the resolved reactions 
jet production from these processes could not be studied in detail in previous fixed target 
photoproduction experiments. At HERA energies however, the results of several authors 0] 
agree in predicting that the photoproduction of jets is dominated by resolved reactions at low 
and medium transverse energies E\ et of the jets, say below 30 GeV. Thus, in the framework of 
this QCD picture, measurement of the jet cross section can be used to obtain information on the 
photon's parton content assuming the parton densities in the proton are known. From e~f deep 
inelastic scattering experiments, studied at e + e~ colliders, the quark content of the photon is 
already relatively well known in the fractional momentum range of the parton 0.007 < x 7 < 1 |J. 
The photoproduction of jets at high energies therefore offers a new tool for the determination of 
the gluon density in the photon. Recent experiments which study jet production in 77 scattering 
are also sensitive to this quantity ||. 

The present analysis studies 2-jet production with more than 7 GeV transverse energy E 3 t et 
per jet in photon-proton scattering. The scattered electron is tagged at small angles so that the 
photon is almost real and the energy of the photon is known from = Et, eam — Et ag - Inclu- 
sive differential cross sections da/dpt and da/drj are derived at the leading order parton level. 
Here pt describes the transverse momenta of the partons and rj their pseudo-rapidities in the 
laboratory system. The unfolding method described below extracts partonic cross sections from 
the measured jet distributions utilizing a Monte Carlo model which besides the LO QCD model 



describes the influence of initial and final state parton showers, multiple parton interactions, 
hadronization and detector effects. 

The results can be directly compared to LO calculations using different parametrizations for 
the parton densities. In next to leading order (NLO) QCD a jet algorithm has to be introduced 
at the parton level as well. The calculations presented so far (?], ||] take a cone algorithm and 
find differences between LO and NLO predictions of the order of only 10% — 30% for values of 
transverse jet energy, jet pseudo-rapidity and other jet parameters used in this paper. 

Besides comparing the inclusive jet cross sections with theoretical predictions, the gluon 
density in the photon can also be derived in a more direct way. The momentum fraction 
x 7 = -E( parton / 7 ) /-E 7 can be fully reconstructed knowing the energies and angles of the jets 
and the energy of the incoming photon. By unfolding the measured x^ ls distribution to the 
parton level using the same Monte Carlo model mentioned above and correcting for the quark 
contribution and the direct photon processes, it is therefore possible to determine the gluon 
density distribution in leading order. 

The paper is structured as follows: after a short description of the detector and of the event 
selection, the Monte Carlo model used for comparisons with the data is studied with respect 
to the energy flow around jets. The model is then used to extract from the jets observed in 
the data inclusive parton cross sections at the leading order QCD level which can be directly 
compared to analytical QCD calculations. It is further used to determine a distribution of the 
fractional momentum x^ from 2-jet events which is interpreted in terms of a) the direct photon 
contributions, b) the resolved contributions with a quark from the photon, and c) the resolved 
contributions with a gluon from the photon. 



2 Detector Description and Selection of 2-jet Events 

A detailed description of the HI apparatus can be found elsewhere Q. The following briefly 
describes the components of the detector relevant for this analysis, which makes use of the 
calorimeters, the luminosity system and the central tracking detector. 



The LAr calorimeter []10|] extends over the polar angular range 4° < 9 < 153° with full 
azimuthal coverage, where 9 is defined with respect to the proton beam direction (+z axis). 
The calorimeter consists of an electromagnetic section with lead absorbers, corresponding to 
a depth of between 20 and 30 radiation lengths, and a hadronic section with steel absorbers. 
The total depth of the LAr calorimeter varies between 4.5 and 8 hadronic interaction lengths. 
The calorimeter is highly segmented in both sections with a total of around 45000 cells. The 
electronic noise per channel is typically in the range 10 to 30 MeV (1 a equivalent energy). 
Test beam measurements of LAr calorimeter modules have demonstrated energy resolutions of 



a(E) /E 0.12/ VE © 0.01 with E in GeV for electrons flu) and a(E)/E « 0.5/V£ © 0.02 
for charged pions pL O]. The hadronic energy scale and resolution have been verified from 
the balance of transverse momentum between hadronic jets and the scattered electron in deep 
inelastic scattering events and are known to a precision of 5% and 10% respectively. 

The calorimeter is surrounded by a superconducting solenoid providing a uniform magnetic 
field of 1.15 T parallel to the beam axis in the tracking region. Charged particle tracks are 
measured in two concentric jet drift chamber modules (CT), covering the polar angular range 
15° < 9 < 165°, and a forward tracking detector (FT), covering the range 7° < 9 < 25°. 

The luminosity system consists of two TICl/TlBr crystal calorimeters having a resolution 
of a(E)/E = 0.1/yE with E in GeV. The electron tagger is located at z = —33 m and the 
photon detector at z = —103 m. The electron tagger accepts electrons with an energy fraction 



between 0.2 and 0.8 with respect to the beam energy and scattering angles below 9' < 5 mrad 
(6' = ir-6). 

The events used in this analysis were taken during the 1993 running period, in which HERA 
collided 26.7 GeV electrons on 820 GeV protons, and correspond to an integrated luminosity of 
290 nb" 1 . They were triggered by a coincidence of the electron tagger and at least one track 
from the central jet chamber trigger. Events were selected, if they fullfilled the following criteria: 



1. The energy deposited in the electron tagger was in the range 8 < Et ag < 20 GeV. The 
cross sections refer to a scaled photon energy of 0.25 < y < 0.7 and a negative squared 
four-momentum of the photon of Q 2 < 0.01 GeV 2 . For the sample used to determine the 
inclusive parton cross sections an additional containment cut for the electron shower was 
applied in order to facilitate the acceptance calculation of the electron tagger. 

2. At least one track in the central tracker with transverse momentum above 0.3 GeV coming 
from the interaction region was required to determine the position of the vertex along the 
beam axis. 

3. The width of the vertex distribution along the beam axis was a = 10 cm. Events were 
accepted in the region of ±3 standard deviations around the nominal vertex position. 



Jet reconstruction was based on purely calorimetric measurements using a cone algorithm 13] in 
a grid of the azimuthal angle ip cel1 and pseudo-rapidity rf el1 which extends from — 3 < rj cel1 < 3. 
The cone radius R = \f Ar/ 2 + Aip 2 in the standard analysis was chosen to be R = 1.0 and 
R = 0.7 was also used for cross checks. Jets were ordered according to the transverse energy in 
the cone. Events were accepted, if 



1. At least 2 jets were found, each with transverse energy above E\ et > 7 GeV. 

2. The jets were contained in the LAr calorimeter < rf et < 2.5. 

3. The rapidity difference between the two most energetic jets was less than |Ar?| < 1.2 in 
order to reject events where the photon spectator (Fig.[l|b) is misidentified as a jet from 
the hard parton-parton scattering process. 



The trigger efficiency was determined to be 94±1% using a monitoring trigger with full efficiency 
for this event selection. A correction of the jet energy scale with respect to the Monte Carlo 
model leads to an additional selection criterion which is described below. The total number of 
2-jet events remaining is 366 without the cut on the shower containment in the electron tagger, 
and 292 for all cuts. 



3 Monte Carlo Generator for QCD Processes 



For the analysis of the data, the PYTHIA 5.6 event generator for photon-proton interactions 14] 
was used together with a generator for quasi-real photons. PYTHIA is based on leading order 
(LO) QCD matrix elements and includes initial and final state parton shower models. The 
strong coupling constant a s was calculated in first order QCD using Aqcd = 200 MeV with 4 
flavours. The renormalization and factorization scales were both set to the transverse momentum 
Pt produced in the parton-parton scattering. Since the QCD calculation used here is divergent 
for processes with small transverse momenta of the partons emerging from the hard interaction 
a lower cut-off has been applied in PYTHIA. This was set to pt > 2 GeV. 



For the proton structure the GRV-LO [15] leading order parton density parametrizations were 
used. The GRV-LO [16] leading order parametrizations were used for the photon structure. The 
latter give a consistent description of the data as will be shown below. Optionally, PYTHIA 
allows for additional interactions within the same event. These are LO QCD processes between 
partons from the photon remnant and partons from the proton remnant. This so called multiple 
interaction option has been explored in proton-antiproton collisions before [O, 17] and the same 



parameters have been used here. For the hadronization process the LUND fragmentation scheme 
was applied (JETSET 111 



The detector response for the generated events was simulated with a detailed simulation 
program and then reconstructed with the same program as used for the data. The generated 
events therefore allow the calculation of correlations between the jets reconstructed in the detec- 
tor and the underlying parton kinematics (see below). These correlations will be used to study 
the parton-parton scattering processes with the jets observed in the data. 



4 Energy Flow in Jet Events 

The precision of the measurement of the transverse jet energy El et and how well the jet energy 
correlates with the parent parton momentum pt are critical matters. This is because the trans- 
verse jet energy distributions fall like (E-J. et )~ n where n ~ 5.5 [Q] so that an imperfect description 
of the energy flow around the jet direction by the Monte Carlo model is a potential source of 
serious error in the conclusions on the parton scattering process. The fact prevented relevant 
conclusions on the photon structure in previous publications [jl], ||, Q and is taken into account 
in this paper for the first time at HERA. 

In Fig.^a the observed transverse energy flow around the jet direction per event is shown 
versus the rapidity distance from the jet axis in a slice of \ip cel1 — (f^ et \ < 1. As an example 
jets were selected with transverse energy 7 < E{ e 1 < 8 GeV collected in a cone of size R = 1 
and pseudo-rapidity between < rf et < 1. The jet profiles are asymmetric, showing a higher 
energy level in the direction of the proton (Ar/ > 0) compared to the photon direction (Ar/ < 0). 
It is interesting to note that the energy flow depends not only on rapidity, but also on the 
momentum fraction x 7 of the parton from the photon side. The parton momentum fraction can 
be reconstructed using the two jets with the highest transverse energy E\ et in the event and 
their pseudo-rapidities rf et together with the energy of the photon E 1 : 

The jet profiles are shown in two bins of x™ s , above and below 0.4. Two effects are observed: 

1. In the photon direction (Ar/ < 0) the low x™ s data show an enhanced energy flow relative 
to the high x^ s data. This can be ascribed to the remnant particles of the photon which 
should be reduced at high x 7 , and absent altogether in the case of direct photon processes 
for which x 7 = 1. 

2. In the proton direction (Ar/ > 0) the energy flow is also enhanced in the low x]^ s distri- 
bution compared with the high x™ s data. This increased energy flow indicates additional 
event activity for events with small parton momenta x 7 , or large energies ~ (1 — x 7 ) of 
the photon spectators. 

The jet core is narrower in the case of the high x^ ls sample compared to the low x™ s sample. 
This effect is connected to the energy flow around the jet and the jet energy interval: since the 
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Figure 2: Observed transverse energy flow versus the rapidity distance from the jet direction 
(integrated over \ip ceU — ip> et \ < 1) for a) data and b) PYTHIA with multiple interactions. The 
jets were selected with transverse energy 7 < E\ et < 8 GeV in a cone of size R = 1 and jet 
rapidities between < r) 3et < 1. The energy flow depends on the momentum fraction x 7 of the 
parton from the photon: shown are large x™ s > 0.4 (open circles) and small x v ™ < 0.4 (filled 
triangles) . 



jets were required to have ~ 7 GeV in the cone, a reduced energy flow around the jet enforces 
more energy in the core. 

In Fig.Qb the transverse energy flow is shown for events of the PYTHIA model including 
multiple interactions. The model also shows an excess of low x^ ls over high x™ s data at positive 
as well as at negative rapidities. 

The observed transverse energy flow in the azimuthal angle |Ay?| around the jet direction 
is shown in Fig.||a,c for two bins of the jet rapidity. The energy flow has been integrated in 
a slice of \rj cel1 — rf et \ < 1 around the jet axis. The transverse energy of the jets was again 
restricted to 7 < E\ et < 8 GeV summed in a cone of size R = 1. The HI data are shown as 
full circles. The observed energy flow outside of the jet cone is much higher for jets at large 
rapidities 2 < rf et < 2.5 (Fig.||c) compared to jets in the central detector region < rf et < 0.5 
(Fig.||a) and is increasing in between 0.5 < rf et < 2 (not shown). Due to the cut on the rapidity 
difference of the two jets |Ary| < 1.2 a part of the second jet is always seen at Aip ~ ir. 

The average transverse energy flow per cone area nR 2 determined outside of the jet cone 
provides a measure of the jet pedestal, often called underlying event [jn]]. This jet pedestal 
is later used to estimate the amount of energy inside the jet cone which is, according to the 
PYTHIA model with multiple interactions, not due to the fragmentation of the hard scattered 
partons. The pedestal energy was determined event by event in the slice \ff el1 — rf et \ < 1 and 
> 1. Regions affected by the second jet were excluded from the measurement. 
The distributions of the jet pedestals corresponding to the jet profiles of Fig.[|a,c are shown in 
FigJ|b,d. 

The PYTHIA model without additional interactions (dashed line) does not give a good 
description of the jet profiles and pedestal distributions. In this model the pedestal energy 
corresponds to QCD radiation and fragmentation effects. The QCD radiation effects are ap- 
proximately taken into account by the parton shower models included in PYTHIA. 

The PYTHIA model with multiple interactions provides an improved description of the jet 
shapes, which is reasonable in the rapidity range < rj :,et < 0.5, but shows deviations from the 
data at large jet rapidities 2 < rf et < 2.5. The differences between data and Monte Carlo model 
are still too large to be neglected. Moreover the sometimes large contributions of energy due 
to multiple interactions in the jet cone give a poor correlation between parton and jet energies 
because this contribution leads to events entering the Ef e > 7 GeV sample which have parton 
transverse momentum below the cutoff used in the Monte Carlo model, namely 2 GeV. 

In order to avoid this problem the transverse energies of the jets in the data were reduced 
by the average pedestal difference between data and PYTHIA without multiple interactions. In 
this way the jets were only corrected for the additional energy flow over and above the initial 
and final state radiation and fragmentation effects. A corresponding pedestal subtraction was 
applied to the jets of the PYTHIA events with multiple interactions. The corrections were 
parametrized in terms of the rapidity of the jets so as to match the effect seen in the data. The 
corrections vary between 0.3 — 2.3 GeV for the data and 0.3 — 1.3 for the events of the PYTHIA 
model with multiple interactions. This subtraction evidently also corrects for the difference in 
the pedestal energies between Monte Carlo and data, provided that this energy difference is the 
same on average inside and outside the cone; which would be the case for multiple interactions. 
As an additional safeguard the transverse energy of the jets is required to be above E\ et > 7 
GeV even after this jet pedestal correction which ensures that the parent parton transverse 
momentum exceeds 2 GeV. In all figures of the following sections these jet energy corrections 
will be applied to data events and Monte Carlo events. 

Overall, the multiple interaction option gives an improved description of the energy flow and 
offers a natural explanation of the observed effects. However, the improved description of the 
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Figure 3: a,c) Observed transverse energy flow versus the azimuthal angle with respect to the 
jet direction (integrated over \rj cel1 — rf et \ < 1) in two rapidity bins: a) < rf et < 0.5 and c) 
2 < rf et < 2.5. The jets have transverse energies 7 < E\ et < 8 GeV in a cone of size R = 1. 
Full circles are HI data. The histograms refer to PYTHIA events with (full line) and without 
(dashed line) multiple interactions. b,d) Distribution of the transverse energy measured outside 
of the jets in the slice \rf el1 — rf et \ < 1, normalized to the area ttR 2 . 



jet shape by this model cannot be considered as a definite proof of the existence of multiple 
interactions. Therefore a model dependence remains in this analysis where the PYTHIA version 
with multiple interactions is used to extract direct information from the observed jet spectra on 
the underlying (LO) parton-parton scattering processes. 



5 Jet-Parton Correlations 



The predicted correlations obtained with PYTHIA between jet and leading order parton quanti- 
ties are shown in Fig.||: a) for azimuthal angles 99; b) for pseudo-rapidities rj; c) for the transverse 
jet energy E\ et and the parton transverse momentum pt, and d) for the reconstructed parton 
momentum fraction from the photon side (Equation ||) with respect to the true x 7 at the parton- 
parton scattering process. In the leading order picture, two hard partons exist per event, and 
there are at least two jets. Therefore for the (p, 77, and pt distributions the jet with the highest 
E\ et was correlated with the parton giving the smallest invariant jet-parton mass, and the second 
jet with the other parton. 

Also, in the case of x™ s (Figjp), the two highest E j t et jets were used. According to the 
nominator of ([!]), different values of x™ s correspond to different 2-jet configurations: small 
values of x™ s require two jets with large rapidities and small transverse jet energies. Here the 
remaining jet pedestal, discussed in Section 4 (dashed line of Fig.||c), raises x v J s relative to the 
true Xry. Events with large x™ s have at least one of the jets at small rapidity or large transverse 



jet energies. At small rapidities the remaining jet pedestal is small relative to the total Ej et 
(dashed line of Fig.||a), and x™ s corresponds in average to the true x 7 . At x 7 ~ 1 transverse 
energy deposited outside of the jets reduce x™ s relative to the true x 7 . The photon energy in the 
denominator of (|l]) was precisely determined from the energy of the scattered electron measured 
in the electron tagger system. The resolution in the logarithm of the reconstructed parton 
momentum fraction is approximately Gaussian in /o<7io(x 7 ) — logio(x'!^ s ) and varies between 
0.22 at true parton momenta around x 7 ~ 0.05 and 0.16 at x 7 ~ 0.5. 

From the study of the jet-parton correlations, the momenta of the hard partons can be 
reconstructed from the measured jet energies. It is then possible to extract partonic cross 
sections and leading order parton densities. The cross sections can be compared directly to 
leading order QCD calculations instead of smearing the theoretical cross section with transverse 
momenta of the partons inside the proton and the photon, QCD radiation effects, fragmentation, 
detector effects, jet formation, and then confronting it with the data. 

In principle such an unfolding procedure seems to be straightforward. The distribution of an 
observable g det measured by the detector is related to the distribution of a partonic observable 
part an i n tegral equation which expresses the convolution of the true distribution with 
all effects between the creation of the hard parton and the measurement process g det {u) = 
J A(u,w)f part (w)dw. This integral equation can be transformed to a matrix equation. Solving 
this matrix equation thus leads directly to the histogram f part (w) and therefore e.g. to the 
partonic cross section. This simple method can produce spurious oscillating components in the 
result due to the limited detector resolution. Therefore the method has to be improved by a 
regularization procedure which reduces the resulting correlations by optimizing both the number 
and position of the bins for the unfolded variable. 



In the analysis presented here an unfolding method [20| has been used to determine the 
differential cross section dajdpt and to derive the distribution of the parton momentum x 7 . 
Because the resolution in the transverse momentum pt is much worse than the resolution in 
rapidity r/, for the extraction of da/drj the transverse momentum pt is unfolded for three different 
bins in r\. The pt integrated results then give the inclusive cross section da/dr}. 
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Figure 4: Correlations between reconstructed jets and leading order parton quantities accord- 
ing to PYTHIA with multiple interactions: a) azimuthal angle (p, b) rapidity 77, c) transverse 
parton momentum pt (to account for the steeply falling distribution the histogram has been 
weighted with the transverse jet energy (E^) 4 ), d) momentum fraction x 7 of the parton from 
the photon side. 



Technically, the unfolding procedure delivers a weight for each Monte Carlo event in terms 
of the true value of the observable chosen for the unfolding which then can be used to reweight 
Monte Carlo distributions of different variables. The comparison of these distributions with the 
data give an important check of the transformation, as described in the Monte Carlo simulation, 
between the partonic distributions and the measurable distributions. An accurate Monte Carlo 
description of this transformation is essential for a reliable unfolding. 

6 Parton Cross Sections 

In Fig.|5] the unfolded single inclusive parton ep cross section da/dpt is shown. The integration 
over the parton rapidity < r/ < 2.5, the initial photon energy and momentum transfer squared 
has already been carried out. The unfolding was done in the pt variable alone. Migration effects 
in pt influence the result in the lowest pt bin most strongly. In order to minimize these effects an 
option in the algorithm was used which allows to constrain the cross section in the bin 4 < pt < 7 
GeV (not shown in the figure) to a reasonable extrapolation of the results shown in the figure. In 
addition, it ensures small bin to bin correlations (less than 30%) and a smooth behaviour of the 
reweighting function. Because the pt distribution is steeply falling, the largest systematic error 
stems from an error in the determination of the calorimetric energy. Variations of the constraint 
on the cross section in the unseen bin and other unfolding parameters result in a 20% error on 
the cross section shown in the figure. Because the result is given in the parton rapidity range 
< r] < 2.5, but the correlation used for unfolding is for all parton rapidities, an additional 
systematic error of 25% is included to account for a possible difference in the correlations of the 
full and the restricted sample. Other systematic error sources are described in Section 7. The 
error bars are computed from the quadratic addition of statistical and systematic errors. The 
cross section, statistical and systematic errors can be found in Table |l[ 

The solid line in Fig.|| is the result of a leading order calculation pl[ ] using the GRV parton 
densities in leading order QCD for the proton and the photon. The factorization and renor- 
malization scales are given by p\. The QCD parameter A was set to 200 MeV. Leading order 
calculations for other photon parametrizations (LAC1: dashed line, LAC3: dotted line |2^]) are 
also included. A NLO QCD calculation has been carried out by the same authors |2l| using 
the GRV higher order parton densities in the photon and the proton, the 2-loop calculation of 
a s , A^g^ = 200 MeV, and a cone size of R = 1. Compared to the LO calculation the NLO 
calculation is essentially larger by an overall factor 1.25. 

As discussed in the previous Section an important check of the success of the transformation 
consists in the comparison of a set of measured observables with the result of the reweighted 
Monte Carlo simulation. This is shown in Fig|]. The distribution of the measured transverse 
jet energy E 3 t et is well described by the Monte Carlo simulation after reweighting (solid line). 
Because the initial set of parameters already gives a good description of the data (dashed line) 
the new weights differ only slightly from one. Also the agreement between the data and the 
model in other variables like rf et for all jets and rf etx + rf et<i for the two highest E\ et jets (Figs. 
^ b,c) is at a satisfactory level. 

The unfolded inclusive parton ep cross section da/dr] is shown in Fig together with the 
theoretical predictions. Statistical and systematic errors are treated in the same way as in Fig.[|. 
The flat shape of the distribution is well reproduced by the LO QCD calculations using the three 
parametrizations of the photon structure function as above. The absolute rate is consistent with 
the calculation using the GRV parametrization. The NLO QCD calculation again is larger by 
an overall factor 1.25 with respect to the LO calculation. 
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Figure 5: Single inclusive partem ep cross section unfolded to the leading order partonic matrix 
element versus transverse momentum pt for < rj < 2.5. The solid line is the partonic cross 
section obtained from a leading order QCD calculation [^TJ using the GRV leading order parton 
densities for the proton and the photon. For the dashed (dotted) line the photon parton densities 
are taken from the LAC1 (LAC3) parametrization. 
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Figure 6: Distributions for the unfolding of pt- Shown is the comparison of data on detector 
level with the Monte Carlo model after (solid line) and before (dashed line) reweighting. The 
parton densities are taken from the GRV parametrizations. a) The E\ et distribution of the jets, 
b) the pseudo-rapidity distribution of the jets, c) the distribution of rf etx + rf et2 for the two 
highest E 3 t et jets. 
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Figure 7: Single inclusive parton ep cross section unfolded to the leading order partonic matrix 
element versus pseudo-rapidity rj for pt > 7 GeV. The solid line is the partonic cross section 
obtained from a leading order QCD calculation [Ell using the GRV leading order parton densities 
for the proton and the photon. For the dashed (dotted) line the photon parton densities are 
taken from the LAC1 (LAC3) parametrization. 



7 Parton Momentum Distributions in the Photon 



In this Section the full 2-jet kinematics is used to check the validity of the leading order QCD 
description of the hard scattering process and to extract information on the subprocesses which 
contribute to the observed jet rate. This will finally allow the determination of parton densities 
in the photon as a function of x 7 . The main goal is to derive, for the first time, the gluon 
momentum distribution a; 7 g(x 7 ) in the photon over a large range in x 7 . 

The correlation between true and reconstructed parton momenta is shown in Fig.^|d. The 
migration effects are large, so the unfolding procedure - described in Section 5 - is used to correct 
from the reconstructed x™ s to the true x 7 distribution. As explained before, the unfolding 
enforces agreement with the distribution (Fig.||a), whereas the agreement of the reweighted 
distributions of other variables simulated by the Monte Carlo model with the actually observed 
data distributions (Fig.|8]b-d) delivers important checks of the transformation: 

1. Fig.||b) the fraction of the parton momentum from the proton side was determined using 
Xp ts = 0.5 (E( ett "e^" + E\ et2 e 1 ^" 12 ^ /E p . For most of the events Xp ts is above 0.01 where 
the quark content of the proton is well established by lepton-nucleon scattering experiments 
and the gluon content is known to a level of 15%. 

2. Fig.||c) the difference in the rapidities of the jets is related to the scattering angle 0* of the 
partons in the CM frame of the parton-parton scattering process: tanhd^ 6 * 1 — rf et2 \/2) = 
cos0*. The distributions of cosG* are predicted by LO QCD for all combinations of the 
interacting partons. This is a basic QCD prediction which has to be fullfilled. 

3. Fig.||d) the differences in the transverse energies of the jets are large compared to the 
experimental resolution and reflect large average transverse momenta of the incident par- 
tons. Initial state radiation effects are important and to a smaller extent large transverse 
energies of the photon remnant. 



The distributions in Fig.||b-d demonstrate that a leading order Monte-Carlo model for the 
hard scattering process is able to give a consistent description of the observed 2-jet events by 
adjusting only the photon structure function. Higher order effects as described by initial state 
parton showers and multiple interactions to describe the enhanced energy flow at large rapidities 
are however essential to describe the data. 

Fig.|9] shows the unfolded x 7 distribution compared to calculations of the PYTHIA generator. 
The full line represents the prediction of the resolved photon processes with a quark on the 
photon side at the hard parton-parton scattering process. The quark distribution in the photon 
as determined by two photon experiments is input in the form of the GRV-LO parametrization 



[16]. The quark induced processes contribute mainly to the central region of the distribution 



corresponding to parton momenta around x 7 = 0.5. 

The prediction of the direct processes are shown as a dashed line above x 7 > 0.77. They 
account for 3/4 of the events in the highest x^ bin. Together with the quark component the size 
of the direct photon calculation, using the GRV-LO parametrization of the parton densities in 
the proton, is consistent with the data at large x 7 . 

At small parton momenta x 7 < 0.2 the contribution of events initiated by a quark from the 
photon as predicted by the given model is clearly below that seen in the data. This suggests 
explaining the additional events at small x 7 by gluons entering the hard process from the photon 
side. 



By subtracting the predicted direct and resolved quark contributions from the data the event 
distribution of the gluon contribution remains. Comparison with the distribution calculated for 
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Figure 8: Distributions of the 2-jet data (circles) are compared with the PYTHIA calculations 
before (dashed line: absolute prediction) and after (full line) the reweighting procedure: a) 
momentum fraction of the parton from the photon x™ s , b) momentum fraction of the parton 
from the proton Xp ts , c) difference in the jet rapidity \rf etl — rf ettl \ which is a measure of the 
scattering angle in the parton-parton system, and d) the difference in the jet transverse energy 
\Ef etl - Ef et2 \. 
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Figure 9: 2-jet event distribution (points) of the true fractional momentum x 7 of the parton 
from the photon. Only the statistical errors are shown. The full line represents the predicted 
contribution of the quark resolved photon processes, the dashed line shows the size of the direct 
photon calculation as obtained by the PYTHIA Monte Carlo. 



processes initiated by a gluon from the photon side, yields correction factors to be applied to 
the modelled distribution. The resulting gluon density is shown in Fig]l(] (see also Table |l]): the 
total errors include statistical and systematic errors where all contributions have been added in 
quadrature. The systematic errors were calculated from the following sources: 

1. Error in knowledge of the calorimeter energy scale of ±5%. 

2. The statistical error in the determination of the hadronic pedestal correction of ±10%. 

3. Uncertainty in the quark density of the photon estimated conservatively as ±30%. 

The systematic errors are dominated at small x 7 values by errors |l| and ^ whereas at high values 
of x 7 error |3] is as large as errors |j] and |2[ An additional error from the luminosity measurement 
of 5% is not included in the FigJTo|. 

To check whether or not a gluon content in the photon is needed to explain the observed 
rate in the data, the following test was carried out using the x v J s distribution: the observed 
data rate in Figj9| was reduced by lowering the energy scale of the calorimeter according to the 
error source |]. The remaining sample was further reduced by subtracting more energy from 
each jet according to error source |2[ The Monte Carlo prediction for processes with a quark 
from the photon however was raised according to the uncertainty given in item || above. With 
the assumption that no gluon exists in the photon the remaining data events were compared 
in a x 2 test to the predicted number of quark and direct events of the Monte Carlo calculatio 
giving a probability below 0.1% for this hypothesis. Thus, a gluon contribution in the photon 
is needed to explain the observed data 

The average parton transverse momentum of the selected events is < pt > 2 = 75 GeV 2 where 
Pt is used as factorization and renormalization scale for the QCD calculation. A change of A 
from 200 MeV to 300 MeV reduces the parton density by 30%. A variation of the gluon density 



in the proton by ±15% [ 23 1 changes the gluon density in the photon by =Fl5%. 



The gluon density in the photon in FigJT(] is compared to the parametrizations of GRV- 



LO[16|, LAC1 and LAC3[22|. A high gluon density at large parton momenta, as suggested 



by the LAC3 parametrization, is clearly excluded. This is consistent with previous observa- 
tions at HERA @, TRISTAN and LEP ||. Above x 7 > 0.08 both the GRV-LO and LAC1 
parametrizations of the gluon distribution are consistent with the data. The strong rise of the 
LAC1 parametrization below x 7 < 0.08 is not supported, while the GRV-LO distribution is 
consistent with the data. In the latter model the gluon distribution of the photon is essentially 
generated by QCD radiation. 



8 Conclusion 

Photoproduction processes with at least two high E\ et jets in the final state were studied with 
the HI detector. The jet quantities as observed in the detector were unfolded to leading order 
parton quantities using jet-parton correlations based on a specific version of the PYTHIA QCD 
generator. Within the present experimental errors the data are - after a correction of the jet en- 
ergy - consistent with this leading order description of the hard scattering process together with 
higher order processes in the initial and final state and multiple interactions, i.e. interactions 
between partons of the photon and proton spectators. 

Single inclusive differential ep parton cross sections were derived as a function of the trans- 
verse momentum pt > 7 GeV and pseudo-rapidity < ij < 2.5 of the partons for photon energies 
corresponding to 0.25 < y < 0.7 and momentum transfer below Q 2 < 0.01 GeV 2 . They agree 




Figure 10: The gluon density of the photon divided by the fine structure constant a = 1/137 
(data: circles) at the scale < pt > 2 = 75 GeV 2 . For comparison the GRV-LO (full) and the 
LAC1 (dashed) and the LAC3 (dotted) parametrizations are shown. 



Pt [GeV] 


8.2 


10.5 


13.9 


19.6 


da/dp t [nb/GeV] 


3.57 


0.90 


0.31 


0.049 


&stat 


0.19 


0.09 


0.05 


0.013 


®syst 


1.52 


0.32 


0.13 


0.018 



rj 


0.42 


1.25 


2.08 


da/dr] [nb] 


5.61 


5.25 


3.50 


&stat 


0.37 


1.02 


0.16 


&syst 


2.10 


2.07 


1.60 





0.059 


0.14 


0.33 


0.59 


0.93 


x^ r g(x 1 )/a 


1.92 


1.19 


0.26 


-0.12 


-0.08 


&stat 


0.87 


0.34 


0.24 


0.15 


0.61 


&syst 


1.68 


0.59 


0.33 


0.33 


0.30 



Table 1: Single inclusive differential parton cross sections (Fig.^,0) and the gluon distribution 
in the photon (Fig,|lO|) together with their statistical and systematic errors. 

well with QCD predictions. A direct study of 2-jet kinematics was used to determine the momen- 
tum fraction of the partons from the photon. At large x 7 it is consistent with the predictions of 
the direct component plus the resolved component initiated by quarks from the photon. At low 
parton momentum fraction the observed jet rate in the data can only be explained by a gluon 
component in the photon. For the first time a leading order gluon distribution was derived down 
to x^ = 0.04. The average scale was < pt > 2 = 75 GeV 2 corresponding to the mean transverse 
momentum squared of the final state partons. Gluon distributions with high density at large x 7 
or steeply rising gluon distributions at small disfavoured. 
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